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G^ ' 

a^ 

0^ ' The static current-current correlation leads to the definitional zero frequency divergence (ZFD) 

in the optical susceptibilities. Previous computations have shown nonequivalent results between two 
(~| ■ gauges (p ■ A and E • r) under the exact same unperturbed wave functions. We reveal that those 

^ ' problems are caused by the improper treatment of the time-dependent gauge phase factor in the 

optical response theory. The gauge phase factor, which is conventionally ignored by the theory, 
is important in solving ZFD and obtaining the equivalent results between these two gauges. The 
Hamiltonians with these two gauges are not necessary equivalent unless the gauge phase factor is 
properly considered in the wavefunctions. Both Su-Shrieffer-Heeger (SSH) and Takayama-Lin-Liu- 
D ' Maki (TLM) models of trans-polyacetylene serve as our illustrative examples to study the linear 

; . I susceptibility x through both current-current and dipole-dipole correlations. Previous improper 

results of the x calculations and distribution functions with both gauges are discussed. The im- 
portance of gauge phase factor to solve the ZFD problem is emphasized based on SSH and TLM 
jrt i models. As a conclusion, the reason why dipole-dipole correlation favors over current-current corre- 

lation in the practical computations is explained. 

^ . PACS numbers: 78.66.Qn, 42. 65. An, 72. 20. Dp, 78.20.Bh 

O 

I. INTRODUCTION 

> . ^^ . . . . 

lO ' The static current-current correlation (JqJo) PI^rI is widely applied in the the optical response theory for many 

decades, both in the definition of linear susceptibility x or conductivity cr^^' , |2] and in the definition of the nonlinear 

optical susceptibilities such as x^"\ where n>2. Mm 

Within the semiclassical theory of radiation, which is also literaturely emphasized by Mahan, |0,H] Bucher and 
^>^ ■ Cotter, B Bloembergen, 0] Shen, and Mukamel, [R the electric field are treated classically and propagation of 
^\ . electromagnetic waves in a medium is governed by Maxwell's equations, that is, the electric field E(r, t) at some 

specific position r and time t can be described as: g] 
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E(r,t)=Eoe*''-'"~'"', (l.f) 



' ^ where Eq is the amplitude, k and lo are the wave vector and frequency. Then, the x in ^ miit volumn v under the 
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static current-current (Jo Jo) correlation is conventionally defined as follows: |H,p| 
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where i7 = — > uji, n^ '^ and m are the electron density and electron mass, eo the dielectric constant, / a unit dyadic, 
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5n,i Kronecker symbol, and 
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y"drdie-*'^'-+''"*(fJo(r,t)Jo(ri,ti)...Jo(r„,t„)), (1.3) 

with T is the time-ordering operator and Jq is the static current operator. H 



If we choose the static dipole-dipole (DD) correlation, the nth-order susceptibiHties will be obtained as follows: 
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ydrdte-"'"-+*^"(fD(r,t)D(ri,ti)...D(r„,i„)), (1.4) 

where D is the static dipole operator. 

It is commonly held that the apparrent zero frequency divergence (ZFD) through the static current experssion 
Eq.(1.2) is only a virtual problem, and the gauge E • r and p • A will give the exactly same results with the same 
unperturbed wave functions. M In other words, Eq.(1.2) and Eq.(|1.4D should reach the same result if one proceeds 



properly. In linear response theory, for the homogeneous and isotropic medium, Martin and Schwinger have shown 
that the ZFD in conductivity a^^' through current expressions could be cancelled by introducing diagmagnetic term 
(A^ term). [|l^ The cancellation of ZFD in the linear conductivity by diagmagnetic term is also discussed by Mahan, 
Haug and Jauho in their famous books, |2|Jl 1[| with a careful consideration on the limitation sequence between k and 
Lu. Ip^ In solid st ate , for either full-filled or empty bands, Aspnes had shown the equivalence of two gauges (from 



Eq.(L2) and Eq.(1.4)) based on the assumption of the cancellation of a ZFD term Eq.(2.6) in his seminal work of 
x") computations. |]l3| These works strengthen the common feelings of the equivalence between two gauges under the 
static dipole or the static current expressions. And the issue of ZFD to be a virtual problem was usually borrowed 
with limited justification. 

Although the 'virtual' property of this ZFD are literally emphasized by a lot of people, strictly speaking, as to 



our knowledge, the solution to this ZFD problem is seldom directly obtained from JqJq correlation [like Eq.(1.2)] 
except Martin and Schwinger's original proof, jl^l As we can see in the calculation of x*^^^ in solid state, the ZFD 
term is conventionally isolated from the convergent term and then is discarded without careful direct check. M] 
Historically, this ZFD problem is of no interest due to the following facts: (i). The gauge transformations seem to 
guarantee the equivalence of two gauges under the same set of wave functions; (ii). In transport theory, 0,M the 
correct imaginary part of Jo Jo correlation (Re[cr'^^-'(ijj)]) still can be obtained and (iii) the ZFD problem (related to 
Im[(T^^^(w)]) usually can be avoided by applying the Kramers-Kronig (KK) relations on the imaginary part of JoJ 







correlation. [pjTl,lJ,nq] Thus the equivalence of two gauges becomes a well-accepted opinion in the community and 



the ZFD is widely considered to be at most a complex technical problem. PJl^] 

Besides those assumptions made in the previous ZFD proof, |0,ttQ-^13] however, in the practical applications, there 
always exists some puzzles from the above common feelings. In a careful study on the gauge choice of two-photon 
transitions in ls-2s transition of hydrogen atom, B Bassani, Forney and Quattropani have found that by directly 
applying the exact same unperturbed wave functions, E • r converges much faster than p • A by using a limitted 
number of discrete intermediate states; Numerical computations also show a 50 % difference of the transition rate 
between those gauges if we only included all discrete intermediate states. Thus, Bassani et al. draw the conclusion 
that E • r is a better gauge. While the study of oscillator strength in a superlattice, by Peeters et ai, turns out that 
E • r and p • A is nonequivalent in the barrier region based on numerical computations. ||17| Thus, on the contrary, 
Peeters et al. concluded that position operator r in the solid state should be redefined and p • A is much better. A 
recent study on zinc-blende semiconductor ||l^ again brings the questions on the equivalence role of E • r and p • A 
in representing the transition matrix of x formula. [|l9| 

Extensive studies on the optical properties of polymers |l|,po|-p3| are based on the tight-binding approximate (TBA) 
models, such as Su-Shrieffer-Heeger (SSH) model |gj| and Takayama-Lin-Liu-Maki (TLM) model [g5[ in weakly 
correlated systems, Hubbard and Pariser-Parr-Pople (PPP) models in strongly correlated systems. Those models 
drastically reduce the complicity of the systems and provide a reasonable way to reveal the actual physical insights 
of many-body systems. In providing the gauge invariance, a U(l) phase transformation has been suggested p3,Eq] in 
those models. Based on the static current formula without diagmagnetic term, there will be ZFD problem in the linear 
conductivity (t^^'{ijj). However, this ZFD problem has not been pointed out clearly and obviously neglected by the 
previous works. po| , p3[ To avoid ZFD in cr^^^ and to obtain a convergent result, Batistic and Bishop suggest to subtract 
the term {[jo, jo]){uj — 0), 122] which is supposed to be the diagmagnetic term directly derived from Hamiltonian. 
Unfortunately, as we will show in this paper that the diagmagnetic term derived from U(l) transformation p3[ can 
not directly return the expected {[jq, jq])(uj — 0) term. Moreover, experimentally observed two-photon absorption 
peakj2^,pq] in the x'^'' spectrum of trans-polyacetylene has raised wide interest for the theoretical explanation. From 
Eq.(1.2) based on TLM models, two-photon cusp was obtained analytically, fit] but it has been criticized by the others 



from the dipole formula approach and other physical concern. p9|-|35| Recently, a quite different analytical form of the 



X^^^ spectrum ||3^ is obtained under DD correlation from that under Jq Jo correlation. ||l| Those above discrepancies 
indeed already cast some doubts on the rooted belief of equivalence between Jq Jo and DD correlations. 

In this paper, we will reexamine the concepts of the gauge transformations and directly show that the gauge 
phase factor, which is not sufficiently emphasized previously and is conventionally ignored under the current-current 
correlation scheme, [nlp| is actually very important in the optical response theory to solve this ZFD difficulty, and 
to recover the equivalent results between the two gauges (p • A and r • E). Therefore, the static current operator 
Jo is no longer suitable for considering the equivalence between two gauges, instead, we should include the induced 
field currents (IFC) which are introduced by the gauge phase factor. To illustrate the effect of gauge phase factor 
rather than the abstract concepts, we only choose the linear response under the one-dimensional (Id) periodic TBA 
models like SSH and TLM models as our examples. Certainly the concepts of the gauge phase factor under the 
specific linear examples also can be expanded into the nonlinear optical response theory and two-dimensional (2d) or 
three-dimensional {3d) cases. 

The paper is organized as follows: In Sec. II, we will reexamine the concept of the gauge transformation and discuss 
the importance of the gauge phase factor in the optical response theory. The problems caused by the ignorance of the 
gauge phase factor are discussed under a general scheme independent of the models. To give an intuitive picture, linear 
optical response under periodic tight-binding-modcls (TBA) such as both SSH and TLM models is investigated in this 
paper, the Hamiltonian under DD correlation (E • r) is discussed (Sec. III. A), and we will study the finear susceptibility 
X^^) through DD correlation in Sec. III. B. The qualitative different results by the conventional polarization operator 
P in optical response are discussed based on periodic models in Sec. III. C. x under current-current correlation is 
discussed in Sec. IV. The SSH Hamiltonian of p • A will be obtained in Sec. IV. A. Before applying the gauge phase 
factor, the ZFD problem of x through Jo Jo correlation is illustrated in Sec. IV. B, previous qualitatively different 
solutions of x^^^ and the practical difficulties in this ZFD problem based on the models are also analyzed (Sec. IV. B). 
After applying the gauge phase factor to the wave functions, a convergent result can be obtained and the ZFD 
problem will be solved (Sec. IV. C). The conditions of equivalence between two gauges are discussed in Sec.V.A and 
the influence of gauge phase factor on the initial distribution function /n(k;) under two gauges are investigated 
(Sec.V.B). The reasons of some previous puzzles on the choice of the gauges are discussed in Sec.V.C. The conclusion 
emphasizing the implications of our work will be described in Sec. VI. 

II. GAUGE PHASE FACTOR IN GAUGE TRANSFORMATION 

Gauge transformation has already been well-understood in the optical response theory. fS^JS^ The equivalence 
between two gauges are built up on the concept of gauge transformation. 
If the electromagnetic field applied, the Schodinger equation is given by: 
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V'(r,t), (2.1) 



where il){v^t) is the exact wave function at space position r and specific time t, m is the particle mass, q is the 
electrical charge, V{y) is the potential, A and (j) as vector and scalor potential correspondingly, under the following 



transformation: 



A^ A' = A + V/(r,i) 

0-0' = 0-|/(r,t), (^-^^ 



where f{r,t) is arbitrary, and A' and (j)' are new vector and new scalor potentials after the transformation Eq.( p.2[ ). 
It could be shown [Q that the form of the Schodinger equation will be exactly the same if the old wave function ijj 
makes the following change into the new exact wave function ■0': 

ijj^i};' ^ e'-^s^'-^*)^ = fair, t)il}, (2.3) 

where gauge phase factor Fg{r,t) is defined as: 

Fg(r,i)^|/(r,t). (2.4) 



The above Eq.(2.2) and Eq.(2.3) are called gauge transformation (or C/(l) transformation pq|). 



Long-wavelength approximation gj^ is used in this paper, that is, the k = in Eq.(l.l), and the electric field E is 
described as E = Eqc"*"*. 

If we consider the following initial scalor and vector potentials under E • r gauge: 

A==0,(/) = -E-r. (2.5) 

After chosing the gauge phase factor Fg as 

F, = '^, (2.6) 

by Eq.(^.2|), we obtain the new vector and new scalor potential under p • A gauge as: 

A' = -,0'-O. (2.7) 

lUJ 



The connection between the old and new wave function is determined by Eq.(2.3) 



Under perturbative schemes to study the optical response, conventionally people use the exact same set of unper- 



turbed wave functions ip'^(r,t) of Hamiltonian Hq (when A = and ijj — in Eq.(2.1)) to serve as our expansion basis 
for both E • r and p • A gauges. [PPP] However, we should point out that the wave functions for both E • r and p • A 
gauges (before and after gauge transformation) should also be restricted by the gauge phase factor Fg from Eq.( p.3[ ), 
therefore two basis sets for both gauges are not the exact same unperturbated wave functions ipni^^ ^); but are different 
by the gauge phase factor Fg. And the Hamiltonian under two gauges(E • r and p • A) are not necessary equivalent 
if they are treated independently and are isolated from the connection between the wave functions under the two 
gauges. Unfortunately, this crucial point has not been clearly illustrated and obviously missed by the perturbated 
schemes works. y,^||§l Especially under current-current correlation scheme, the gauge phase factor's contribution is 
obviously ignored and A'^{t) term is considered as no any physical meanings. [p9[ Thus the current-current correlation 



is conventionally reduced into the JqJq formula such as Eq.(1.2), and the equivalence between current-current and 
dipole-dipole correlations is usually considered as Jq Jo ^^^d DD correlations under the exact same basis of unperturbed 
wave functions. jl|,|||,|l|,|l| 

An elegant review by Langhoff, Epstein and Karplus covered the topics of time-dependent perturbative theory, BO] 
they have sharply pointed out that the time-dependent phase in wave functions is very essential and the improper 
treatment of time-dependent phase will cause secular divergence in time-dependent perturbations. In field theory, 
it is also well-understood that the improper treatment of the phase factor will cause divergence. 0| Since the gauge 



phase factor Eq.(2.6) is obviously time-dependent, neglecting this phase factor will cause the ZFD in the susceptibility 



computations, as the examples we will show in the following sections. 

III. LINEAR RESPONSE BY DIPOLE-DIPOLE CORRELATION 

Having appreciated the importance of the gauge phase factor, we choose the following single electron periodic 
models - SSH (or Hiickel) model and TLM model as our examples provided the following reasons: (i). Those periodic 
models are widely applied in the polymer theory in the 80s and early 90s; Remarkable results have been obtained; |4l| ] 
(ii). The optical susceptibilities obtained from those models can be analytically solvable and be compared with the 
previous results, [g^,^ (iii). In both the SSH and the TLM models, Peirels instability ^^ leads to the semiconductor 
property of two band structure - with valence band full-filled and conduction band empty. It is very obvious from 
the physical point of view, as the freqency of the electrical field goes (reaches the static electric field), the linear 



conductivity ct'^^ and the linear susceptibility x will not be reduced to Drude formula as in the metals, |14|l5[| and 
will not cause ZFD problem. In this section, we will first discuss x^^-' and cr^^-' of infinite chains (with number of {CH) 
unit N goes to infinity) under the DD correlation. 

A. SSH Hamiltonian 

Based on periodic TEA, The SSH Hamiltonian [|4| is given by: 

. to + i-iy^ 

Ls '- 



HsSH = - / , 



{Cl,^A.s + ClCi+us), (3.1) 



where to is the transfer integral between the nearest-neighbor sites, A is the gap parameter and Cj ^{Ci^g) cre- 
ates(annihilates) an tt electron at site / with spin s. In continuum limitation, above SSH model will give the TLM 
model. 12^ For the SSH model, each site is occupied by one electron. 

If we want to include the electron-photon interaction E • r directly from the polarization operator P, where 



p = J2RicIci 



(3.2) 



and 



Ri = /a+(-l)'u 



(3.3) 



is the site I position with the lattice constant a and dimerized constant u, p4| we will face the problem of ill-definition 
of P in periodic systems. [pJ3| , |l7|j2^ ] To solve this problem, we should consider the imposed periodic condition of the 
position operator r. fl^Jl^l Expressing position operator r under the Bloch states |n, k >= u„_k(r)e''''", where w„,k(r) 
is the periodic function under the translation of lattice vector, [Q we will be able to satisfy the periodic condition of 
r as follows: 



and 



where v is unit cell volume. 



Tnk.n'k' = «(5„,„'Vk5(k - k') + n.n,n'{^)S{k - k'), 



f^n,«'(k) = - / <k(r)VkM„',k(r)rfr, 



(3.4) 
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We change Hamiltonian Eq.(3.1) into the momentum space by applying the following consecutive transformations: 
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where 



s{k) = \J [2tocos{ka)f + [l\sin{ka)Y 



(3.9) 



and Ri^ and Ri^ are odd and even position defined by Eq.( |3.3[ ). a^git) and a^ g{t) are the excitations of electrons in 
the conduction band and valence band with momentum k and spin s. 

If we choose the spinor description -0^ s(i)= ('^fc s(^)' '^fc s(^))' ^^^ Hamiltonian including E • r in momentum space 
is described by: 



HssH{k,t) ^Ho + Het, 



(3.10) 



where 






and 



i?E.r = -£> • Eoe'- 



By Eq.(3.4), the dipole operator D could be obtained as foUows: 
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where 



p(k) 



72(fc) 



U, 



(3.14) 



is the coefficient related to the interband transition between the conduction and valence bands in a unit cell 2a and 
the second term in Eq.( 3.13| ) is related to the intraband transition Jl^ e is the electric charge and a are the Pauli 
matrixes, u is dimerized constant related to the lattice distortion. 1241 



B. Linear response through E ■ r 



For the hnear susceptibility, Xssh^^-''^'^) '-^^ ^^ obtained from Eq.(L4) and Eq.( |3.13 ): 
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(3.15) 



where the Green function G{k,uj) is the fourier transformation of G{t — t') = —i{T^jj{t)t(j'' {t')) as follows: 

^,, s UJ+UJk(T3 

G[k,uj) = ^ TT"' 



(3.16) 



with LUk = e{k)/h and e = 0+. 

By Eq.(|3l5|), we have XslHi^)= x'ssh(-^^^)- 



(1) _ eH2toa) n 



(1 — rjSx^ydx 



[(l-<52a;2)(a;2-l)]*a;2(a;2_z2)' 



(3.17) 



where x = hujk/^, z = huj/{2A), S = A/(2io) and relative distortion r/ = {2u)Ja. Eq.( 3.17 ) can be numerically 
integrated if one change x — > a; + ie in considering the life-time of the state. |2 9^32(1 For polyacetylene, by choosing 
to = 2.^eV, A = 0.9ey, a = 1.221, u = 0.041 and e - 0.03, (2|-|3|] we have 5 = 0.18 and i] = 0.07. The value 
of Ixsiifl with or without 77 contribution are plotted in Fig.l. As we can clearly see from the graph, the relative 
distortion 77's contribution is very small (about 1%). We can see the 'unklapp enhancement' peak at z = 1 comparing 
with the peak z = 1/6. Those results are also discussed in the previous works. pS] 




,(1) 



FIG. 1. \XssHi'^)\ with z = hLu/{2A), for e = 0.03, S = 0.18 and for rj = 0.07 (solid line) or for r; = (dashed line). 



If the continuum limitation is applied, that is, S —f 0^ , rj —> 0^ , e —f 0^ and 2toa —f hvp, the above integral 
Eq.( |3.17 ) approaches the linear optical susceptibility Xtlm(^) under the TLM model |2^] as follows: 
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where 
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cosh^ (z) iiT 

zy/z'^ - 1 2zVz2 - 1 



(z^ < 1), 



(z2 > 1). 



(3.19) 



The conductivity a{uj) given by —iuj ■ x*-^-*- The i?e[cr'^'(w)] is the exact same as previous results. [pO|j2 
The calculated Xtlm ^^'^ absolute value of Xtlm ^^^ shown in Fig. 2 and Fig. 3. 
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FIG. 2. The real part (solid line) and the imaginary part (dashed line) of XTLAii^) with z = hii}/{2A) 
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FIG. 3. 



IXTLAfC"^)! with z: 



hbj/{2/\). 



The above computations indeed are based on perturbative scheme with ■0]^ g{t) and il'k,s{t) as unperturbed creation 



and annihilation operators under Haniiltonian Hq. [ Eq.( |3.11 )] From the above figures and expressions, it is very 
obvious that the DD correla tion approach (or E • r) will not have ZFD in x^'^'' and a'^^K The straightforward compu- 
tations easily show that Eg. ( 3. 17 ) and Eq.( 3.18| ) obey the KK relation. Those results are certainly reasonable under 
the physical picture. 



C. The qualitative different results through the polarization operator P 



If the above dipole operator D is replaced by the polarization operator P defined by Eq.(3.2), which are extensively 
used in the models, we will obtain a different result from Eq.(3.18). As pointed in the Hteratures, p^ , p^ , p7| , p3[ 
the polarization operator P is sensitive to boundary conditions and singular in the thermodynamic limit. We will 
demonstrate this issue based on P in SSH model. Although the original point of position Ri can be arbitrary, we can 
pick up the chain region I from 1 to A^ units to do a simple test. By Eq.(3.2), Eq.(3.6)-Eq.( |3.9D , we obtain the unit 
polarization P""**(fc) from the total polarization _P*°*°'(fc) in the momentum space: 



Uk) = lim 
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total 



k.s 



(3.20) 



We find that P^™*{k) contains no intraband transition like 2nd term in Eq.( ^.15| ). Applying a similar computation 
as DD correlation, we can obtain linear susceptibility Xssh based on P""'*(fc): 
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(3.21) 



where F(7r/2, Vl — 5'^) is the complete elliptic integrals of the first kind. Obviously, the 1st term in the above Eq.(B.21) 
is a constant and it is logarithm singular if J — > 0^! This result obeys the previous conclusions obtained in solid states, 
|0,n6J17[23] showing that the polarization operator's failure in optical susceptibility's study. We choose the same 
parameter as in Sec. III. B to plot spectrum of Ixsiff I' ^^^ comparison curves plotted in Fig. 4 showing that the results 
computed from D and P are actually different both in the spectrum and in the magnitude. We find another peak at 
z = 1/5 corresponding to the transition from the bottom of valence band to the top of conduction band is greater 
than the peak at z = 1. The 'unklapp enhancement' disappears although the same position of resonant peak still can 
be predictable based on P. Those differences show that the electron transition from all the other sites can not be 
neglectable especially in the periodic models. And conclusions obtained in x on the dependence of gap parameter 



A are qualitatively different for an infinity chain between D and P. |45| Therefore there are some nonconsistencies 
between the polarization operator P and the dipole operator D at least in the linear response theory. 
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FIG. 4. \XssH ('^)l under polarization operator P (solid line) versus \XssHi^)\ under dipole operator D (dashed line) with 
z = huj/{2A) and i] = 0. There is also a magnitude difference of 10^ between D and P besides the qualitative differences in 
the spectrum. 

IV. LINEAR RESPONSE THROUGH CURRENT-CURRENT CORRELATION 
A. SSH Hamiltonian in the vector potential form 



T he t ight-binding Hamiltonian under p • A form should be invariant under gauge transformation (Eq.( |2.2D and 
Eq.(^.3[)). p3,p6[ If we change the phase of the one-particle wave function in tight-binding approximation; 



Ci{r)=e^<^Cs{r). 



(4.1) 



we must modify the kinetic energy term according to the unperturbed Hamiltonian Hq in the Wannier functions basis 
as follows: 



Ho^Y. ^(r-r')(Ct(r)a(r')+Ct(r')a(r)) 
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+ C:t(r')e'''r''''-^^"^/''C:(r) 



(4.3) 



where t(r — r') is the hopping from position r to r', C](r) creates an electron at site r with spin s, q is the particle 
charge. The above transformation is also known as Peierls substitution, p^ The above form Eq.( |4.3| ) has some kind 
of general meanings in TEA and can be frequently seen in the theoretical works, p^ 



If the function f{v,t) is arbitrary in Eq.(2_^), it is easy to verify that the above TEA Hamiltonian (43) is an 
invariant if the local phase change according to Eq.( |4.l|) or Eq.( ^.3|) : 

e{v.t)^^f{v,t). (4.4) 

As a specific example, the SSH Hamiltonian with the vector potential A should be as follows (we change A into A 
since it is Id case): 1431 



Hssh{A) = -Y, 



l,s 



'to + {-1)4 



(C";+i.e-''='^(^'^-^'+i)/'"'C" 



l,s 



+ C'l,e''''^<-^'-^'+'^/'''C'i+i^s), 



(4.5) 



B. result without gauge phase factor 



Under the assumption of the gauge phase factor (4.4) as no physical meanings, p] we ignore the phase 9 and treat 
the creation and annihilation operators {C; } and {C'l} the same as unperturbed creation and annihilation operators 
{C/} and {Ci} ||,|6| defined in Eg Q. 

The above p • A Hamiltonian Eq. ([4.5|) can be expanded in the powers of the external vector potential A and obtain 
the following: 



Hssh{A) = Ho- JoA - -JiA^ + OiA^), 



(4.6) 



where Hq is given by Eq.(3.1) 



^o = -E^ 



to + i-iy 



[a-2(-l)'wl 



(C;+l,sQ,s " Ci,,Ci+i^s) 



(4.7) 



and 



■^--sd) 



Ls 



to + i-iYj 



[a-2{~iyu] 



icl,A.s + clA+i,s), 



'1+ 

The current operator J is obtained from the following equation: 

i 



J 



P.H 



From Eq.(3.2) and Eq.(4.6), we obtain the current operator under the SSH Hamiltonian as follows: 

JssH ^ Jo + JiA, 



(4.8) 



(4.9) 



(4.10) 



where Jo and Ji are defined by Eq.(4.7) and Eq.(|4.8|) 



Similar to the computations in Sec. Ill, we transform the Hamiltonian Eq.(4.6) and the current operators Eq.(4.7) 
and Eq.(4.8) into the momentum space. By applying Eq.(|3.6|) and Eq.(3.7), we can obtain the following: 



M^) - y E [Ao(fc)^I_,(t)CT3^fc,.(i) +So(fc)^i.,,(t)'Ti^fc,.(i) 

k.s 



and 



Ji 



(fc) = (^)2 ^ \A,{k)i;ljt)cj^i,uAt) + B,{k)i7^At)a,i^uAt) 



(4.11) 



(4.12) 



with Ao(fc), i3o(fc), ^i(fc)and i?i(fc) defined as follows: 

{2tof{l-5'^)sin{2ka) 



Ao(fc) 



2e(fc) 



2inA 



B,ik)^-^+,eik) 



and 



Ai{k)^{l + ri'')e{k) 



e{k) ' 



Bi(fc) 



{2to)^r){l-S'^)sini2ka) 



(4.13) 



(4.14) 
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r], S and e{k) in the above equations are defined the same as in the Sec. III. 

Applying Eq.( |l.2|) and Eq.(1.3), and considering the diagmagnetic current Ji{k)A obtained from the above SSH 
Hamiltonian Eq.(4.6), we have the linear susceptibility under JqJo correlation as follows: 



x'siff(-'^i''^i) 



>^'L!o(-'^i'^i) , x')J{-^i,u;i) 



(4.15) 



with 



x1o7o("^l'^l) 



E 



(y) 



Trl Ao{k)a3Gik, u)Ao{k)a3G{k, to - wi) 

+ Ao{k)a3G{k, uj)Bo{k)aiG{k, u - uji) 
+ BQ{k)aiG{k, Lu)Ao{k)a3G{k, uj - uji) 

+ Bo{k)criG{k, uo)Bo{k)(JiG{k, uj - ui) 



duj 



(4.16) 



ifi~^.-.)--^qfT.jy\ 



Ai (k)(Tj,G(k, ijj ~ ui) 



i3i(fc)aiG(fc,o.-c^i)l^, 



(4.17) 



where Green function G{k,uj) is defined by Eq.( 3.16|) . Following the straightforward computations, we can obtain the 
following: 



xi:U-) 



-2i 



e'^{2toa) 



(1 — rjSx ) dx 



and 



x'5fH = 2z 



Trr 7i [(l-(52a;2)(a;2-l)]*(a;2-z2)' 
e2(2ioa) /■* [(1 + 772)^23.2 _ 2ri5]dx 



■kTi 



[(l-<52a;2)(3.2_l)]^ 



(4.18) 



(4.19) 



Eq.( 4.19 ) is similar to the first term in Eq.( 3.21 ) and is a constant independent of z. Eg. ( [4. 19] ) obviously is not the 
term ([jo, jo]) ( a; = 0) suggested by Bati stic and Bishop. |2|] 
By Eq.(|4.15|), Eq.([4.18D and Eq.(p^), we obtain: 



,(1) e2(2ioa)J p 

X ssh{^) - 2^A2z2 



(1 — rjSx^ydx 



7 [{l + r]^)d^x^ -2r]6]dx 
[(l-(S2a;2)(a;2_l)]i(a;2_z2) _/^ [{1 ~ S^x^){x^ - 1)]^2 



(4.20) 



where x,A, z,(5 and r/ are all the same defined as in Eq.(3.17) 



If the continuum limitation (same as in Sec. Ill) is applied, we will find that the contribution (Eq.(4.19)) from the 
diagmagnetic term ( Ji term) disappears. We will obtain the following susceptibility under TLM models: 



XTLAli^) 



e^fivp 



27rA2 



/(^), 



(4.21) 



where f(z) is defined by Eq.( |3.19[ ). 

We plot lx'55/fl, X TL M ^^"i Ix'tla/I ^^ Fig-5, Fig. 6 and Fig. 7 with the same parameter as in Sec. III. From the 
analytical form Eq.(4.21) and the Figures, we find that the first term in Eq.( 3.1S|) disappears in this JoJq formula, 
however, this important feature has not been reported in the previous work. [E0D23] 
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FIG. 5. |x'ssh('^)I through JqJo correlation (soUd Une) versus |Xss£r('^)l through DD correlation (dashed line) with 
z = huj/{2A), for e — 0.03, f] — 0.07 and S — 0.18. |x'g5^('^)| obviously shows the ZFD if the gauge phase factor is not 
considered. 
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FIG. 6. The real part (solid line) and the imaginary part (dashed line) of X TLMi'-^) with z = huj/{2/S). It shows the ZFD 
in real part when the gauge phase factor is not considered. 




FIG. 7. {yi tLm{^)\ through Jo Jo correlatoin (solid line) compared with |XTi,i\/('^)l through DD correlation or J J correlation 
through gauge phase factor (dashed line) with z = huj/{2A). |x'ti,m('^)I shows ZFD. 



U nder TLM mo dels, t he missing diagmagnetic term can be understood directly for m the Lagrangi an. |44| ] Obviously, 
Eq.( [4.20| ) and Eq.( [l.2l| ) strongly diverge when z ^ for the real part in both Eq.( fl.2C| ) and Eq.( |4.2l| ). The above 
results are certainly wrong since it does not follow the KK relation besides this ZFD problem. Careful comparisons 
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are made between Eq.(4.2C), Eq.(4.21) and Eq.(3.17), Eq.( 3.18| ), showing that the correct imaginary parts obtained 
from the Jo Jo correlation Eq.(1.2) and Eq.(1.3) are still maitained. It also shows that based on Jo Jo correlation, the 
absorption part (related to imaginary part) will still be correct. This conclusion is not novel to the peoples working on 
the transport problem, where the Kubo formula based on Jq Jo correlation is already applied as a common knowledge. 
B Peoples usually resort KK relation to avoid real part ZED difficulty. [M 

But how to understand those difficulties (such as ZED and violation of KK relation in Eq.( 4.20| ) and Eq.( [4.2l| )) 
becomes a task. The long wavelength approximation used in the above examples already eliminate the possible ZED 
caused by the limitation sequence between k and w. [|l2| Someone maybe argue to solve the ZED problem by including 
a diagmagnetic term, such as the effective mass m* and electron density hq assumption be made in this solid state 
problem. |lq,nq] In fact, those assumptions are awkward because the parameters m* and uq in diagmagnetic term 
cannot be predictable in the models. Eor example, in the above SSH and TLM models, m* and no could be arbitrary. 
pl| In the TLM model, you can not includ e diag magnetic term directly from Lagrangian. El] Eor the SSH model, the 
diagmagnetic term can be included in Eq.(llC) from the TBA Hamiltonian Eq.(L6), but you still cannot solve this 
ZED problem, as we clearly see from Eig.5 and Eq.( 4.20 ). Moreover, based on diagmagnetic terms to cancel ZED, 
the assumption on the property of medium should be made in the proof. P, |lO| , [ll| 

The conventional way to treat this ZED problem is to separate the divergent terms from the convergent term 
and directly throw it away with some possible physical explanations. |13|| Eortunately, in linear response, neglect 
of diagmagnetic term does not cause so many troubles since the imaginary part of Jo Jo correlation is still correct. 
P, [llHl4[| It is not a big surprise that people have not taken this ZED problem seriously although the doubts are always 
existing in those two gauges especially in the models. ||9|, p^ , |38|| 



C. solving ZFD by gauge phase factor 



As we demonstrate above, the diagmagnetic term directly obtained from the Eq.(4.5) can not solve the ZED 



problem. This ZED problem, which is a conceptual problem as we have already pointed out in the Sec. II, is caused 
by the conventional careless treatment of the gauge phase factor in the optical response theory. In this part, we will 
demonstrate the same result of SSH and TLM models as DD correlation after considering the contribution of gauge 
phase factor in J J correlation. 



As we discussed in Sec. II, the new creation operator C; and annihilation operator C'l in Eq.(4.5), should be differred 
by the gauge phase factor from the unperturbed creation operator Cj and annihilation operators Ci. EoUowing 
Eq.( ^.3|) , we obtain the following after the local phase factor is considered: |46|] 

(4.22) 
(4.23) 
(4.24) 
(4.25) 
ij--(fc) = ^^Y,Boik)i^Ut)a,iJ;kAt)- (4-26) 



C'l =p«^«'AC'; 

Thus by the above relation, the Hamiltonian in momentum space should make the following change: 

H{k) ^ H{n), 

where 

, eA 
n 

Eor the SSH Hamiltonian, we have the following new Hamiltonian: 

H"''"{k) = iJ^'"'"(fc) + Hl'^'"{k)A + 0{A^), 



where iJ^'="'(fc) defined as the Eq.(3.11) and 



The new current operator J"'^"'(/c) could be obtained from the commutator equation [D{k),H^'^^{k)]/{ih) as the 
following: 



J«'=«'(fc) = J^'™(fc) + J["''"{k)A + 0(^2), 



(4.27) 



where J^'^^ik) is exact the same as Eq.(4.11) and 
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where 



k,s 

Blik) 



(4.28) 



B5'™(fc) = -Ao(fc) 



2ioA 
e2(fc) 



(4.29) 



where Ao(fc) and Bo{k) are defined in Eq.( [4.13| ). 

Thus, we find that after considering the gauge factor, the new current operator J^'^'^ik) is the same as the static 
current Jo{k), but J"'^'"(fc) is different from the static current Ji{k), we call those current differences between J"'^"'(/c) 
and static currents Jo(fc), Ji{k), etc., as induced field currents (IFC) since they are introduced by the gauge field. 

Through the evolution operator in interaction picture, H it is easy to derive the formula of the Xjj as follows: 



X 



(1) 



Xj.e..j-..™(-a;i,u;i) - — 






(4.30) 



where 



.(1) 

3S" 



^ir^K 



,{-iOi,U)i) 



I.. ■/ — ( 



Tr 



Ao{k)a3G{k, uj)Bo{k)aiG{k, w - wi) 
Bo{k)aiG{k, iu)Bo{k)aiG{k, uj-ljA^, 



(4.31) 



Trj Ar"'(fc)a3G(fc, 



U) — LOi) 



B^-^{k)a,G{k,uj-ujA^, 



(4.32) 



Eq.( [4.3lD gives exact t he sa me result as Eq.( [4.18D , which is computed through JqJq correlation. While the contri- 
bution from jp^™ in Eq.(4.32) can be obtained as: 



x'liL(c.) = 2z5'^2toa) 



Jl" 



nh 



(1 — -qdx^Ydx 



[(l-,52x2)(a;2_i)]ia;2 



(4.33) 



The above term (4.33) actually is {[io,Jo]){^ = 0) suggested by Batistic and Bisho p. [p2[ 

By Eq.( |4.1^ ) and Eq.( [l:.33| ), Eq.( |4.30| ) returns the exact the same result as Eq.( |3.17D and Eq.( |3.18D , which are 
computed through DD correlation. 



V. DISCUSSIONS 
A. equivalence condition betAveen two gauges 

The above examples based on Jq Jo and DD correlations show the importance of the gauge phase factor Fg in the 
optical response theory. The equivalence of the two gauges (E • rand p • A) should not be based on the exact same 
sets of wave functions, but should be different by the gauge phase factor. This crucial point has never been clearly 
pointed out previously. The conventional equivalence between the static current-current (Jo Jo) correlation and the 
static dipole-dipole [DD) correlation is not maintained, instead, we should consider the induced field currents (IFCs) 
which is introduced by the gauge phase factor. For example, in the periodic models, we should do the phase shift 
from k ^ K for the basis wave functions (or through creation and annihilation operators in models) to obtain the new 
current operator J'"^™ for linear case. Thus, the equivalence between current-current and dipole-dipole correlations 
should be understood as the situation when all IFCs are included. This paper's conclusion that the time-dependent 
gauge phase factor Fg can not be ignored under the perturbative scheme is also consistent with Langhoff et aZ.'s 
results. Eotl 
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B. Initial distribution function in two gauges 

In optical response, people are much more interested in the population of the states than the phase factor. Recently, 
femtochemistry experiments [^ are able to reveal the phase factor's effect by the vibrational modes of nucleis other 
than the population effect (or distribution function) of states. To some extend, the phase factor's influence on the 
optical response begin to raise people's interest. Although our treatment of the optical susceptibilities is still under 
the frame of non-vibration nucleis, the gauge phase factor's effect in the theory still can be seen from the above 
computations. 

Another important feature caused by gauge phase factor Fg is the influence on the initial distribution function /n(k;) 
between two gauges. If all particles are in the ground state ip'i of unperturbed Hamiltonian _ffo when the electric field 
E is applied at time t — 0, under the perturbative scheme, we can use the set of unperturbed wave functions {"0"} 
with the initial distribution function 

/g(k) = 1 and /„(k) =0 for all theother(n7^.g); (5.1) 



Previously, we use the exact the same initial distribution function such like Eq.(5.1) for both gauges. 

As we show in the previous sections, the exact wave functions between two gauges should be different by a time- 
dependent gauge phase factor. Those conclusions are also kept under the perturbative scheme. Thus, the initial 



distribution function in two different gauges are not necessary the exact same as Eq.(5^) and should be carefully 
considered on the choice of the basis sets. 

Specifically for E • r, we should choose unperturbed wave functions as the basis to avoid ZFD directly, |Q then the 
initial distribution function could be Eq.(5T) since the initial ground state is ip^. But for p • A, since the new ground 



state should be e^^^ij/l according to Eq.( |2.3| ), there are two ways to set the initial distribution function: (i)If we have 
already considered th e ga uge phase factor in our new basis set, that is, we use {e^^s^Oj ^^ q^j. ^jasis, the distribution 



function is still as Eq.(5^); (ii)If we still use unperturbed wave function {V'^} as o ur basis, we should project the initial 
wave function e^^^ipg on the basis set {■0^} instead of directly applying Eq. (p^) as our initial distribution function. 
The latter way is much more complicated since the initial set of distribution function will be time-dependent. In the 
previous example we show, we use the first way under the local phase approximation. [^6| Both ways to treat initial 
distribution function by p • A are complicated, in this meaning, the gauge E • r is much better. This conclusion has 
actually already been tested in the practical applications. p-[7|Js| Jll||l3| , p9| -p6[ 

C. some previous puzzles 

In understanding those nonequivalent results between the two gauges, some literatures are emphasized by the others. 
Q However, those explanations lacks the direct proof and the effect of the gauge phase factor is ignored. The above 
two subsections provide the qualitative reasons to understand those nonequivalent puzzles between the two gauges 
in the numerical computations. P,[13^,p6| Based on the same set of unperturbed wave functions, |9| the computations 



of Bassani et al. shows that the E • r converges much faster than p • A reveal that the distribution Eq.(5.1) is good 
for E • r but not for p • A case, which is certainly reasonable according to our discussions in Sec.V.B. The x*-^-* 
computation is a special case, ||l^,^ since some terms could be if the symmetry is applied. |^^ Unfortunately, the 
nonequivalent results between Jq Jo correlation and DD correlation are magnified in the x computations. As an 
example, the spectrum of x^'^-' of polyacetylene shows a different spectrum between two gauges. [0,|6| Those difference 
could be understood qualitatively through the gauge phase factor. 

VI. CONCLUSIONS 



From the computations and discussions above, we concluded that Eq.(1.2) based on the static current is improper 
and leads ZFD problem, thus, IFCs generated from the gauge phase factor should be included to solve this difhculty. 
Generally speaking, the Hamiltonians under two gauges are not necessary equivalent unless the gauge phase factor is 
properly considered through the wave functions. Because the choice problems in the initial distribution function and 
the basis sets under two gauges (Sec.V.B), it will be very tedious to do the perturbative computations based on the 
p • A than E • r. If the careful computation is made based on the concept of gauge phase factor, both gauges will 
lead to the equivalent result. Although our computations are chiefly based on a Id periodic model, it could be easily 
seen that the chief conclusions of this paper can be expanded in the 2(i, 3d and other systems based on the general 
illustrations in Sec. II. 
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